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1. Introduction 

There exist several definitions of coherent states (CS) [l]-[3]. They lead to the same result for 
the harmonic oscillator potential and usually to different results for other physical systems. 
Nevertheless, a careful analysis shows that they possesses common properties which can be 
taken to define CS for the very general quantum system [3]. It happens that such a definition 
is very ambiguous. Therefore one can choose between different possible systems the one which 
has a desirable property. For instance, one can look for a state, the wave function of which can 
be expressed in terms of known special functions or even in terms of elementary functions only. 
This simphfies considerably the study of these CS and their use in other applications. In this 
way one was able to construct different systems of CS for the well-known soliton potentials [4] . 
It was also shown that different CS may be related between them with a symmetry operator 
[4] or with the Laplace transform [5]. 

The most essential progress in studying CS has been achieved for systems possessing 
symmetries [2]. If it is not possible to associate a symmetry group with a quantum 
system the problem becomes much more complicated. For instance, after applying a 
nontrivial supersymmetric (or equivalently Darboux) transformation to a Hamiltonian allowing 
a symmetry, the symmetry is usually lost or it is transformed to a nonlinear symmetry [6] . This 
makes impossible to apply group theoretical methods for studying CS. In this respect it was 
suggested to apply the same transformation to the known CS of the initial quantum system and 
treat states thus obtained as coherent states for the transformed system [7, 8] . This approach 
proved to be useful for studying CS for supersymmetrically transformed Harmonic oscillator 
[8] and time-independent singular oscillator [9] potentials. In particular, using correspondence 
between classical and quantum systems, which can be realized just by the technique of CS (see 
e.g. [2]), one was able to give an interpretation of the supersymmetry transformation in terms 
of classical notions [6] . Here we continue this study at the level of the time-dependent singular 
oscillator. 

The time-dependent singular oscillator Hamiltonian 

ho ^ -dl + u;\t)x^ + gx-^ (1) 

plays an essential role in different physical applications between which we would like to mention 
interesting results in molecular physics [10], in optics [11] and in mathematical physics [12]. 
Recently a model of a two-ion trap has been proposed (see [13] where a good literature review 
is also given) based on this Hamiltonian. Exact solvabihty of the Schrodinger equation with 
the Hamiltonian (1) may be, in particular, related with the fact that the symmetry operators 
realize a representation of the su{l, 1) algebra and therefore it can be solved by the method 
of separation of variables [14]. Moreover, for this algebra two systems of CS are commonly 
known [13]. The one may be obtained with the help of the SU{1, 1) group translation operator 
(they are known as Perelomov CS, see e.g. [2]) and the other are eigenstates of the annihilation 
operator [15] (Barut-Girardelo CS [16]). Here we will study the supersymmetric transformation 
of either system. 

The paper is organized as follows. In the next section to fix the notations we give a review 
of some properties of the singular oscillator Hamiltonian and its CS. In particular, basing 
on the Barut-Girardello CS we construct a holomorphic representation of state functions and 
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operators; for a concrete example of a free particle moving with a fixed value of the angular 
momentum equal to two wc compare the time stability of Barut-Girardello and Perelomov CS. 
In section 3 we apply the time-dependent Darboux transformation both to Barut-Girardello 
and to Perelomov CS and study some properties of states thus obtained. Conclusions are drown 
in the last section. 



2. Singular oscillator Hamiltonian 

In this section we summarize briefly what is known about this Hamiltonian and its coherent 
states (see e.g. [6],[9]-[13],[15],[16]) we need further. We would like to stress that although 
almost everything presented here is already published, but from one side we are adopting 
here a different approach with respect to the one usually used and from the other side we 
are bringing together a number of facts spread in numerous literature. In contrast to most of 
papers using different modifications of the method of quantum invariants [17] for solving the 
Schrodinger equation we are using the method of separation of variables [14] . The advantage 
of this approach is that we do not need to be restricted by any quantum mechanical picture. 
Instead, we consider the Schrodinger equation as a second order parabolic differential equation, 
which gives us the possibility to get "nonphysical" solutions we need for applying the Darboux 
algorithm. 



2.1. Solutions of the Schrodinger equation 

Generators of the SU{1, 1) symmetry group, being symmetry operators of the Schrodinger 
equation with the Hamiltonian ho, in coordinate representation 

,2 



A;_ = 2 [a^ - s^gx'^] k+ = 2 [{a+Y - e^gx-''\ , (2) 

ko^l {k-k+ - k+k-) (3) 
are expressed in terms of the harmonic oscillator creation and annihilation operators 

a = edx — ^sx a"*" = —ed^ + ^tx . 

The dot over a symbol means the derivative with respect to time. Parameters e and £ are 
solutions to the equation of motion for the classical harmonic oscillator 

e + Au;^{t)e = 0. 

In particular, when e is complex, e will denote its complex conjugate and they are such 
that is — ee = \i. Casimir operator C is expressed in terms of the parameter g: C — 
I {k^k^ + k^k^) — /cq = 3/16 — (7/4. This gives us the relation between g and the representation 
parameter k: C = k{l — k), g = 3/A + Ak{k — 1). 

Solutions ipn{x,t), n = 0, 1, . . . of the Schrodinger equation belonging to the Hilbert space 
of functions square integrable on semiaxis form a basis of the discrete series irreducible 
representation T^{g) of the group SU{1,1) k — \ -\- -\- 46. The vacuum vector 
'4>Q{x,t) is selected from the set of solutions of the Schrodinger equation by the equations 
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koipolxjt) = kipo{x,t), k_'ipQ{x,t) = 0. All other basis vectors may be obtained by acting with 
the creation operator k+: 

The action of the operators (2)- (3) on this basis is given by 

kotjjnix, t) = {k + n) ^„(x, t) k^tjjn{x, t) = -c^ilJn±i{x, t) (4) 

To solve the Schrodinger equation simultaneously with the eigenvalue equation for ko we 
are using the method of separation of variables. From coordinates (x, t) we go to curvilinear 
coordinates {u,v) defined as: u = X'y~^^'^, v = t were 7 = ee. The choice 

^ = e>^"'P{v)Q{u) (5) 

guarantees the separation of variables. The eigenvalue equation for ko is reduced to the first 
order ordinary equation for P which can readily be integrated to give P = 7"^/^ [e/e)^. So, 
replacing the action of k^ on a solution by the multiplication on the separation constant A one 
obtains from the Schrodinger equation the following second order differential equation: 

Qyy-[\y^ + 9y-^-2X]Q = Q y=\u. (6) 

This equation can be reduced to the equation for the Laguerre polynomials if A = n + fc, 
n = 0, 1, . . ., which just corresponds to the well-known coordinate representation of the basis 
ipn{x, t). If A = —k — norA = /c — n — l,n = 0, l...it can give rise to Laguerre polynomials 
also but of course this does not correspond to square integrable solutions of the Schrodinger 
equation. Finally these solutions have the form 

The lower sign and the choice a — 2k — l correspond to the discrete spectrum eigenf unctions. 

The upper sign permitting a to be both positive and negative corresponds to solutions which do 
not belong to the Hilbert space. For a = 2k — 1 they are such that koipnix, t) = —{k + n)iljn{x, t) 
and for a = 1 — 2k they satisfy the equation koipn{x,t) = {k — n — l)ipn{x,t). The value 
Non — 2^'^~^{n\)^r~^{n + 2k) guarantees the normahzation of the square integrable solutions 
to the unity. For the non-normalizable solutions Nq^ does not play any role. 

For applying the Darboux algorithm we need nodeless solutions. The zeros of the Laguerre 
polynomials are well-known [18]. In the negative semiaxis, a; < 0, L"{x) has only one node 
provided (a + l)„ = {a + l){a + 2) . . . (a + n) < 0. It is clear that since we put k = i + ^Vl + 46, 
a = 2k — 1 = ^y/l + 46 > meaning that in this case any ipn is nodeless and it is such that 
'0~^(cxo) ~ but '^"^(0) = 00. For a negative a the result depends on the interval where a 
falls. If —2m — 1 > a > —2m — 2 we have n = 0, 2, . . . , 2m and for —2m > a > —2m — 1 more 
values of n are possible: n = 0, 2, . . . , 2m, 2m + 1, 2m + 2, . . . where m = 0, 1, . . .. Moreover, 
for a < —3/2 the function ip~^{x) is square integrable and belongs to the domain od definition 
of ho when it is considered as an operator in the Hilbert space. 
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2.2. Barut-Girardelo coherent states 

The states ipfj_{x,t) known as Barut-Girardelo CS are defined as eigenstates of tfie operator k- 

k^'il)^{x,t) ^ H'il)^{x,t) iieC. (8) 

To solve this equation simultaneously with the Schrodinger equation we are using the method 
of separation of variables once again. In coordinates {u = x/e,v = t} the Schrodinger equation 
separates if = exp{jeeu'^)P{v)Q{u). The function P{v) is defined by a first order ordinary 
equation which is easily integrated to give 
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P = s-h-^^ 

The Schrodinger equation is reduced to the second order equation for Q: 

Q"{u)-{gu-^ + fi^)Q{u)=0 (9) 

which after the change of the dependent variable Q = u^I becomes the modified Bessel 
equation. The condition for ilj^{x,t) to belong to the domain of definition of ho selects us 
only one solution to this equation: / = /2fc_i(/iM). (We are using the standard notations for the 
Bessel functions [18] .) To calculate the normalization integral we are making use of the tables 
[19]. So, the normahzed solution is 

Expanding the Bessel function in a Taylor series one can get their expansion in terms of 
Laguerre polynomials and in such a way they are expressed through the basis functions (7) 
with coefficients depending only on even powers of fj,. Therefore, it is more convenient to 
change the complex variable /i in favor of X — 2/1^ which yields 



M^, t) = A^oA V a„A>„(x, t) a„ = ( /|")/^ (11) 
^ ^yn\^{n + 2k 

with 

A^oA := (V'oIV'a) = (AA)^-^ V-i (2 |A|) (2k) . (12) 
As usual the states i/^xix, t) are not orthogonal to each other 

(V^A' I ^a) = ^ ^ r ■ 

[/2.-1 (2 |A|) 72.-1 (2|A'|)]^ 

Since the basis '0„(x, t) is complete in the Hilbert space one can calculate the measure p(A) 
which realizes the resolution of the identity over the coherent states V'a [9] 

j M{^^Jx\dp{X) = 1 p(A) = ^K2k-i{2\X\)hk-i{2\X\)dXdX. (13) 

All integrals over the variable A are extended to the whole complex plan. 
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Now one can construct a holomorphic representation of the vectors and operators [9] . Any 
ip{x,t) — J2'^=o'^n'ipn{x,t) from the Hilbert space of square integrable functions defined on the 
positive semiaxis can be written in coherent state representation '0'^(A): 

oo 

r (A) := {ij-^{x, t) I ij{x, t)) = Nox J2 ^nCnX"" = Noxi^ (A) A e C . (14) 



n=0 



The holomorphic function V'(A) = Yl'^=o '^nCn^"' can be associated with the usual square 
integrable function ip{x,t) given by its Fourier coefficients c„ over the basis (7). 

Using the complex conjugate form of the resolution of the identity (13) one can define an 
inner product {ipa (A) | ipb (A)) in the space of the functions ■0a,6(A) holomorphic in the complex 
plane 

{1/Ja\lpb)= I |iVoA|'^a(A)^„(A)rfp(A):=(^„(A)|^,(A)). (15) 



This means that the integration in the space of holomorphic functions should be carried out 
with the measure dp{X) — |A^oa|^'^p(A), so that 

(A) I i,, (A)) = ji,^ (A) ^, (A) dp (A) . (16) 

To distinguish this inner product from the one in the space L^(0, oo) we indicate the integration 
variable inside the brackets. The space of holomorphic functions il^{X) such that 



/ 



\^{X)Ydp{\) < oo 

equipped with the inner product (16) becomes a Hilbert space. 

The orthonormal basis il^nix^t) in this representation looks hke as follows: 

i^n (A) = a„A" (V'n (A) I Vn' (A)) = 5nn' ■ 

From here we find the Dirac-dclta function 

oo oo 

5 (A, A') = 5] V'n (A) (A') = «n(AA> = F {2k) {XX')'-'%,., (2 VaT) 

n=0 n=0 

which takes off the integration in this space: 



^(A)= J 6{X,X')ij{X')dp{X'). 

The CS in the holomorphic representation 

^A'(A) = (Ay)^-'=/2fe-i(2VAA^) 
and the SU{1, 1) generators 

ko^X^ + k k+^X k^^X^ + 2k^. (17) 
dX dX^ dX 
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are time-independent. The mean values of the operator ko in the coherent state 



and 



(2|A|) 



{^|;, l^ol V'a) ^k' + |Ar + |A| (19) 



may be useful in the following. 



2.3. Perelomov coherent states 



Perelomov CS ijjz{x,t) are obtained by acting on the ground state function with the group 
translation operator. Therefore their Fourier coefficients over the basis ipn{x, t) are independent 
on t and coincide with the ones for the time-independent Hamiltonian 

oo 

xIj,{x, t) = No, J2 M^, < 1 (20) 

n=0 



Using the generating function for the Laguerre polynomials one gets their explicit expression 
^.(..*)=2i-»r-^(2*).-V-^(i^)'exp(-|-i±| + !0) (22) 



where C = z- . 

^ £ 

We notice that since the series expansion (20) of the time-dependent CS (22) in terms of 
the time-dependent basis {ipn{x,t)} is exactly the same as the similar expansion for the time- 
independent Hamiltonian, all properties of such states known for the stationary case take place 
for the non-stationary one. In particular, they realize the resolution of the identity operator 
and, hence, one can map any square integrable on the positive semiaxis function to a function 
holomorphic in the unit disc getting in such a way a holomorphic representation. We will not 
go in details for this case since they are very well-known [2] . 



2.4. Discussion 

Since the parameters z and A labelling the CS fall into different domains, it is not 
straightforward to make a comparison between them. One of the possibilities could be to 
choose these parameters such that the mean value of an operator is the same for both states. 
As a particular example we take uj = (free particle) and g = 2 (orbital quantum number 
1=1). In this case e = -^[t + i) and we have chosen real values for z and A such that the mean 
value of ko coincides with its value in the first excited state (7) which is /c + 1. In this case 
z = {2k + 1)~^/^ and A ~ 1.021. With these values of z and A we plotted the squared modulus 
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Figure 1. Comparison between proba- 
bility distribution for Barut-Girardelo 
and Perelomov CS for time moments 
t = 0: 1 and la; t = 0.5: 2 and 2a; 
t = 1: 3 and 3a; t = 2: 4 and 4a. 
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Figure 2. Comparison between proba- 
bility distribution for Darboux trans- 
formed Barut-Girardelo and Perelo- 
mov CS for time moments f = 0: 1 
and la; t — 0.5: 2 and 2a; t — 1: 3 and 
3a; t ^2: 4 and 4a. 



of il)z{x,t) (22) and il)\{x,t) given by (10) at /i = a/A/2 in fig. 1. It is clearly seen from the 
figure that at the initial time moment the Barut-Girardello CS are much more localized then 
Perelomov CS but the latter are much more stable in time. At times grater then 2 they are 
already very spread whereas the Perelomov CS maintain almost the same localization as at the 
initial time moment. In the next section we shall show that the Darboux transformation keeps 
the different time behavior of two types of states practically unchanged. 



3. Darboux transformation of coherent states 



Time-dependent Darboux transformation [20] can create real potential differences only if the 
Schrodinger equation has at least one solution u = u{x,t) satisfying the reality condition 
(logu/w)^^^ = 0. Solutions (5) satisfy this condition for any real function Q. So, taking 
different real solutions of equation (9) one can, in general, obtain a two-parameter family 
of exactly solvable partners for Hq but only the functions u = ipn{.x,t) (7) produce potential 
differences expressed in terms of elementary functions. Between square integrable solutions only 
ipo{x, t) is nodeless. Unfortunately, it produces a new potential of the same kind as Vq only with 
a different value of g (shape invariance at the time-dependent level). So, to get essentially new 
potentials we have to keep in (7) only upper sign. For a < —3/2 any nodeless function (7) gives 
rise to a new potential. For simplicity we will consider here only the case a > when all square 
integrable solutions for the transformed Hamiltonian hi = —d^ + Vi{x,t) can be obtained by 
acting on corresponding solutions of the initial equation with the transformation operator 

L = Li{t)[d^-u^{x,t)/u{x,t)] Li{t) =exp[2 J dtlm{lnu)^^]. (23) 

The potential Vi{x,t) is expressed in terms of the same function u: 



Vi{x,t) = Vo{x,t) + A{x,t) A{x,t) = -[ln|M(a;,t)|^] 
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Taking one of the function (7) with n — m (which is supposed to be fixed from now on) one 
gets the potential difference A{x, t) — Am{x, t): 

^ ~ 87 ■ 

The function Li{t) is determined by (23) up to a constant which we fix to simphfy subsequent 
formulas so that Li = ^/2rf and the explicit expression for the transformation operator (23) is 

L=v^[a.---^- — -^^]. (25) 

Normalized to unity solutions for the Hamiltonian hi are: = Ni^Lipn, — {n-\-2k-\-'m)~^ . 

Operator L and its formally adjoint L"*" factorize the symmetry operator feg- L~^L = 
ko + k + m. The opposite superposition LL^ is a symmetry operator for the transformed 
Schrodinger equation. If we denote po = LL^ — k — m then the functions (pn are eigenfunctions 
of po and the spectrum of po is identical to the spectrum of ko. Using k± and transformation 
operators L and L"*" one can construct other symmetry operators for the transformed equation: 
p± — Lk±L'^. They are ladder operators for the functions </?„: 

P±(p„{x,t) = {NinN^n±i))~^C^ipn±l{x,t) . 

We notice now that this Hamiltonian gives us an example of the Schrodinger equation, symmetry 
operators of which do not close a Lie algebra but satisfy a polynomial algebra 

\po,P±] = ±P± 

[p-,p+] = 2[k{l -k) + 2po{k + m) + 2pl]{po + k + m) . 
Similar algebra was previously obtained for the time- independent singular oscillator in [6]. 

3.1. Transformation of Barut-Girardelo coherent states 

To obtain CS for the Hamiltonian hi we act with L given in (25) on CS ipx. (f\ = NixLipx. The 
factor Nix being calculated from the formula = ('^aI^o + k + m\ipx) and (18) guarantees 
the normalization of the states (px to unity. Their series expansion in terms of the basis {(fn} 
can be found by acting with the same operator on the series (11): 

00 

cpx^ N^bny^Pn bn^ an{n + 2k + m)^2k + m)-^ (26) 

n=0 

where N ^ {2k + m)^NixNox- 

The states (px = ipx{x,t) thus obtained may be interpreted as coherent states if they admit 
the resolution of the identity operator 

J MMdpiX)^!. (27) 
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Now we proceed to find the measure p(A). We will look for the function p(A) depending only 
on the absolute value of A, |A| = y/x: dp — \h {x) dxdcj) ; the function h{x) is to be determined. 
Therefore, it is convenient to use polar coordinates in the complex plan of the variable A, 
A = 6xp . After being integrated over the variable (j) equation (27) yields 



n\V{n + 2k) 

Prom here it follows that if equation 

7r(n + 2/c + m 



(28) 



r(n+ l)r(n + 2fc) 



POO 

/ a;"W2-i/4 j-1 ^ V^) ho (x) dx = 1 (29) 

Jo ° 2 



is satisfied, than (27) will also take place because of the completeness of the system {(fin}- H 
now we rewrite (29) as 

P / N , r(n + l)r(n + 2fc) 

/ (x) dx = ^—^ r — ^ 30 

Jo ^ ' (n + 2A; + p) ^ ' 



where 



^{x)^\NoxN,x\^h{x) (31) 



we recognize in it a problem of moments on semiaxis (see e.g. [23]). To solve this problem we 
are using the following integral [24] 



^0 



x^f {x) = r (n + 1) r (n + 2k) (32) 



where f{x) = 2x'^ 2i^2fe-i(2\/x). It is not difficult to get the expression for $(x) in terms of 



r>oo 

$(x)=x-+2^-W y-"^-^f{y)dy. (33) 

J X 

Indeed, first we notice that since /c > | we have x$ (x) when x — 0. Therefore, the 
integration in (30) by parts under the condition (32) yields just the right hand side of (30) 
meaning that equations (33) and (31) define the measure p(A). 

The resolution of identity (27) permits us to construct a holomorphic representation for the 
transformed Hamiltonian hi. The Fourier coefficients {c„} of a function (p{x,t) over the basis 
{ipn{x^t)} gives us the same function in CS representation: tp'^ (A) = (v^xlv^) ~ {^) ^-i^d the 
function tp (A) = Y1^=q bnCn.X"' is the holomorphic representative of (px{x, t). Now one can define 
a new inner product 

((^i(A)|(^2(A)) = j{cpi\<pj){cpj\<p2)dp{X)^ j \Nfipi (A) cp2 (A)dp(A) 

in the space of holomorphic functions, which gives us a holomorphic representation of states 
and operators different from that discussed in section 2.2. 
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It is not difficult to see that after the Darboux transformation any basis function ipn (A) goes 
to V5„(A) = {n + 2k + m)^{2k + m)~^ipnW- Therefore, if we want for the functions ipn (A) and 
ifn (A) to be related by the Darboux transformation, we have to put 

ifniX) = {n + 2k + m)-^L{X)ipniX) V^„(A) = {n + 2k + m)-5L+(A)(/?„(A) . 

This gives us the holomorphic representation of the Darboux transformation operators 

L(A) = (m + 2k)-^ [ko{\) - /c - m] L+(A) = (m + 2k)^ . 



3.2. Transformation of Perelomov coherent states 



Once again we act with the transformation operator L given in (25) but now on ipz{x, t) (22) to 
get the Darboux transformed Perelomov CS, Lpz{x,t) = Ni^LtpziXji). Normalization constant 
is easily calculated using the equation {Lipz\Lipz) = {■ipz\L~^Lipz) and the factorization property 



of the transformation operators, N^^ 
basis {f n{x, t)} is 



m + 2A;(l — |zp). Their Fourier series in terms of the 



(x,t) = Ar^^6„z>„(x,t) 



(34) 



n=0 



AT, 



NozNiz{2k + m)2 6„ = a„(n + 2k + m)^{2k + m)"2 . 

Here also the series (34) is exactly the same which has previously been obtained for the 
time-independent oscillator [6]. We have already found [6] the measure, which realizes the 
resolution of the identity operator in terms of (pz, constructed new Holomorphic representation 
for the operators and states, found the Kahler potential and symplectic 2-form meaning that we 
obtained a classical mechanics, which being quantized a la Berezin gives us back the holomorphic 
representation of the quantum system. This procedure can be considered as the one giving rise 
to a classical counterpart of the Darboux transformation valid both for time-dependent and 
time-independent cases. 



3.3. Discussion 
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Figure 3. Comparison between proba- 
bility distribution for Barut-Girardelo 
CS before (solid line) and after (dashed 
line) Darboux transformation. 



Figure 4. Comparison between prob- 
ability distribution for Perelomov CS 
before (solid line) and after (dashed 
line) Darboux transformation. 
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We would like to point out that if the Schrodinger equation has a symmetry algebra, 
this property is usually lost after the Darboux transformation. Nevertheless, if it has ladder 
operators, the transformed equation may also have them. In such a case it is possible to look 
for eigenstates of the annihilation operator and call them coherent states. For the case of the 
time-independent harmonic oscillator potential this approach was realized in [21]. Since ladder 
operators have now two derivative orders more with respect to the initial ladder operators, the 
differential equation they satisfy is rather complicated [22] which makes difficult studying such 
states. Our approach has an advantage that the transformation operator (25) is a simple first 
order differential operator. So, it is very easy to operate with it. Moreover, in such a way one 
can get different systems of CS if they are available for the initial Hamiltonian. Usually different 
systems of CS exhibit different properties [11, 15]. We conjecture that Darboux transformation 
approximately preserves different behavior of different CS. To support this conjecture we plotted 
the transformed Barut-Girardelo CS together with the transformed Perelomov CS on Fig. 2. 
Prom the first sight the difference between Figs. 1 and 2 is practically invisible. To show 
it better we plotted Barut-Girardelo CS (solid fine) together with their transformed version 
(dashed line) on Fig. 3. Fig. 4 shows the Perelomov CS before (solid line) and after (dashed 
line) the Darboux transformation. It is clearly seen from these figures that for both cases the 
Darboux transformation results mainly in a displacement of the curve while its shape is very 
little affected. 

4. Conclusion 

We have shown that acting with the Darboux transformation operator on the known CS of 
the time dependent singular oscillator gives us the states with the similar properties. Thus, 
the ones obtained from Barut-Girardelo CS may be called Barut-Girardelo-like CS while the 
others, which are produced using Perelomov CS, may be called Perelomov-like CS. Each system 
of CS admits a resolution of the identity operator which makes it possible to construct different 
holomorphic representations. A particular example of a free particle in p-state {£ = 1) shows 
that Barut-Girardelo-like CS are well localized at the initial time moment while Perelomov-like 
CS are more stable with time evolution. Such a behavior is a reflection of the similar behavior 
of corresponding states before the transformation. Therefore we hope that the new systems of 
CS may be useful in similar applications where the known systems have proven to be helpful. 
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